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Abstract. We consider the dynamics of a model introduced recently by Bialas, Burda 
and Johnston. At equilibrium the model exhibits a transition between a fluid and a 
condensed phase. For long evolution times the dynamics of condensation possesses a 
scaling regime that we study by analytical and numerical means. We determine the 
scaling form of the occupation number probabilities. The behaviour of the two-time 
correlations of the energy demonstrates that aging takes place in the condensed phase, 
while it does not in the fluid phase. 
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Introduction. Recently a mimbcr of studies have been devoted to the dynamics 
of the Backgammon model, a simple stochastic model which exhibits some of the 
features of glassy systems such as slow dynamics, non-stationary properties of two- 
times correlations, violation of the fluctuation-dissipation ratio, etc. [1-8]. The model 
introduced in [9] is a simple generalisation of one of the models defined in [3] (model 
B), itself closely related to the Backgammon model. In contrast with the latter -or 
with model B- it exhibits, at finite temperature, a phase transition between a fluid and 
a condensed phase. The aim of this paper is to study the dynamics of condensation 
in this model, hereafter referred to as model B'. This study is motivated by the fact 
that, while for model B the nonequilibrium properties such as slow dynamics [3] or 
aging [5] only occur at zero temperature, here they appear in a whole phase, where 
the system condenses. In this work we focus our interest on the scaling behaviour 
of the occupation probabilities and just attempt a short description of the two-time 
correlations. The dynamics of this model may also serve as a source of inspiration for 
the understanding of the dynamics of the Bose-Einstein condensation, for which little 
is known. Finally the dynamics of the original branched polymer model introduced 
by the authors of ref. [9] may have an interest in its own. 

Definitions. Consider a system of N particles distributed amongst M boxes. We 
denote by Ni the number of particles contained in box number i (i — 1, . . . , M), with 
J2i = The energy of a given configuration C = {Ni,N2, . . . , Nm} of the system 
is defined as the sum of the energies of individual boxes £{C) = The 
Boltzmann factor associated to E{k) is denoted by pk- The partition function of the 
system reads ^: 

^M,jv = ^ . . . ^ Piv. . . .Piv„5( ^ AT,, N)=I — ^ [P{z)f. (1) 

JVi Nm ■' ^'^^^ 

The right hand side expression is obtained by using the integral representation (5(m, n) 
= § dz z'^~'^ /2'Kiz for the constraint. Equation (1) shows that the eqinlibrium prop- 
erties of the system depend only on the set {pk}, or, equivalently, on its generating 
function P{z) = J^kPkz''- 

In the present work we restrict our study to the class of models for which E(k) 
behaves logarithmically^ for large k. The function P{z) has therefore a finite radius 
of convergence Zc = 1, with a singularity of the form {zc — z)^~^. For definiteness we 
study model B' defined by taking E{k) = ln(l + fc), hence Pfc = (1 + k)~^ leading to 
a Dirichlet series for P{z) [9]. In parallel we will recall below the properties of model 
B [3] defined by taking E{k) = -5k,o, hence pk = e'^^"'" and P{z) = e^ -|- z/{l - z). 

Equilibrium properties. The equilibrium properties of the models follow simply from 
the previous definitions. In the thermodynamic limit {M,N — > oo), the density p = 
N/M being fixed, the method of steepest descent can be applied to the integral (1). 
The saddle-point equation reads dP{z)/dz = pP{z)/z. The saddle-point value Zs, 

^ Note the difference between the statistics used in eq. (1) and that used in the 
definition of the Backgammon model [1, 7, 8]. In contrast to the latter, here the 
particles are not identified by a label. In this sense they are indistinguishable [3]. 
^ Adding to E{k) a linear term in k plays no role because of the constraint. 
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which is by definition the thcrniodynamical fugacity of the model, is thus related to 
temperature and density. The free energy per box reads —f3f = In P{zs) — plnzg. 
When Zs increases from to Zc, df{zs)/dzs = —ln{zs) is positive, hence f{zs) is 
monotonous. One also finds that in this range, p increases monotonously from to 
Pc = ZcP'{zc)/P{zc)- While Pc is infinite for model B, it is finite in the case of model 
B' as long as /3 > 2 and reads pc — ({[3— l)/({f3) — 1 where ( is the Riemann function. 
This fundamental difference between the two models is a consequence of the behaviour 
of the Pk at infinity and is at the origin of the possible existence of condensation in 
model B'. Indeed when pc is finite, f{p) reaches its maximum at f{pc) = lnP{zc). 
Therefore as long as p < pc the system is "fiuid". When p > pc a, condensed phase 
appears [9]. Thus model B has only a fluid phase for T > 0. 

These two phases are characterized by different forms of the occupation probabilities, 
defined as follows. The probability that a generic box, say box number 1, contains k 
particles is defined as fk = Prob{A''i = fc}, i.e., fk represents the fraction of boxes 
containing k particles. The same definition holds out of equilibrium. The conservation 
of the number of boxes and of the number of particles imposes that J2k fk = ^ and 
Sfe ^fk = P- F^om the definition above, one gets 

Ni Nm \ ^ ) 

In the thermodynamic limit, using again the steepest descent method, one obtains 

k 

fk=Pk-^^, {p<Pc) (3) 
in the fluid phase. In the condensed phase, one has 

fk =^y ip> Pc) (4) 

which is the same as eq. (3) with Zg = 1. The normalisation condition /j. is fulfilled 
by both equations. However the conservation of the number of particles, ^j^. kfk = p, 
holds only in the fluid phase, while it is violated in the condensed phase where this 
sum is equal to Pc- The M{p — pc) missing particles sit in a single box [9]. 

An analogous situation occurs for model B at T = 0. Eq. (3) gives fk = pe^^''-"{l — 
Zs)'^z^-'^, with p{l - Zs) = 1- fo. When T ^ 0, pc ^ 0, z^ ^ 1, hence /o 1. 
Again, in order to restore the conservation of particles, all the particles have to be in 
a single box. 

Definition of the dynamics. The rules defining the dynamics of the models follow 
naturally from their static definitions. These rules were given for model B in [3] . In this 
work we used both the Metropolis rule, more convenient for Monte Carlo simulations 
and the heat bath rule, leading to simpler dynamical equations. Let us first describe 
the former one. At every time step 6t = 1/M two boxes are chosen at random, a 
departure box d, containing k particles, chosen amongst the non-empty boxes, and an 
arrival box a, containing I particles. Note the difference with the Backgammon model, 
where the departure box is defined by choosing a particle at random (see footnote 1). 
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The transfer of one of the particles from box d to box a is accepted with a probabihty 

. f. Pk-iPi+i^ 
mm 1, — 

V Pk Pi 

In the heat bath case, once a particle is drawn, it is put into one of the boxes with a 
probabihty proportional to the equilibrium probability of the resulting configuration. 
Thus this move is accepted with a probability ( J2h^o fi^'^) ^ ■ '^^^ correspond- 
ing dynamical equation for the occupation probabilities is the master equation of a 
random walk for A^^i, the number of particles in the generic box number 1: 

= Mfe+l/fe+l + Afe-l/fc-l(l — Sko) — {/J'ki^ — Sko) + ^k)fk, (5) 

where ~ 1, {k > 1) is the hopping rate to the left, corresponding to Ni = k ^ 
Ni = k — 1, and 

1^1=0 J I Pk 

is the hopping rate to the right corresponding to N\ = k ^ = k + 1. The factor 

1 — 5ko accounts for the fact that one cannot select an empty box as a departure box 
nor can Ni be negative, i.e. A_i = /Uq = 0. In other terms a partially absorbing 
barrier is present at site fc = 0. This random walk is biased, to the right or to the 
left according to whether its velocity — /i^ is positive or negative, respectively. It 
is easy to check that eq. (5) fulfills both conservations of boxes and particles. 

In the stationary state {fk = 0) one recovers the equilibrium results given above. 
The detailed balance condition yields fk+i/fk — A^, the two possible solutions of 
which are precisely those given in eqs. (3) and (4) above, from which one gets 

1 - /o _ ( Zs if p < pc 



Model B yields 



^k - _a - (»j 

l-./o + /oe ' 

Whatever the value of (3 the walk is biased to the left. Therefore, intuitively, no 
condensation is expected in this model, except at zero temperature. In this case, since 
Aj; = 1 if /c > 0, and Aq = 0, the system performs a symmetric random walk with a 
totally absorbing barrier at the origin [3]. Hence for t ^ oo, /o ^ 1, i.e. all boxes 
become empty. Nevertheless, in order to fulfill the conservation of particles, one box 
has to contain all the particles, as already explained above. T = appears thus as a 
critical point. Let us point out that for model B, the Metropolis algorithm and the 
heath bath rule lead to the same dynamical equation for the fk- 

The case of model B' is richer. A simple analysis, and a numerical check, show that 
if /3 is large enough, for small k the bias is to the left, while it is to the right for large 
k. Thus one intuitively expects condensation in this model, for a whole range of values 
of /3. 

Condensation in the scaling regime. A numerical integration of eq. (5), and Monte 
Carlo simulations for /J = 4 and p = 2 > Pc(4) = .110, give some insight in the 
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phenomenology of the dynamics in the condcnsc;(l phase. Three regimes are observed. 
First a transient one, with a rapid reorganisation of the particles in the boxes, leading 
to a situation with a fluid part for small k and the appearance of a condensate, i.e. 
a group of boxes containing a large fraction of the total number of particles. This 
regime is followed by the scaling regime, our main interest in this work, where the 
evolution in time of the condensate is self-similar. The number of boxes containing 
the condensate decreases (although it remains large in this regime). 




Figure 1. Scaling in the dynamics of condensation for model B'. kfkVi, obtained by 
numerical integration of eq. (5), is plotted versus k/\/i for 10 different times varying 
from 100 to 1000. (/3 = 4 and p = 2.) 

Finally at very long times, the condensate reduces to a single box according to a 
non universal process where finite size effects should now be taken into account. This 
sequence of three regimes takes place in a similar fashion in the case of model B. It is 
also reminiscent of the dynamics of coarsening systems [10]. 

In order to describe the scaling regime for we set 

^* (l+ewfeH ) if k<uo/e 



fk^{ pi^r ' , (9) 

e'^g{u){l + 0{s)) if k > uo/e 

where the small scale e{t) is to be determined, and u = ek is the scaling variable, uq 
fixes the separation between condensed and fluid phases and corresponds intuitively 
to the position of the "dip" clearly visible on figure 1. 

The normalisation conditions J2 fk = ^ and J2 ^fk = P lead respectively to 

uo/s oo 

E -^^-^ + / 9iu)du = 0, (10) 



and 

P- Pc 



/ ug{u)du. (11) 
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It is easy to check that, once the hniit t ^ oc' is taken, with uq fixed (hence uq/s ^ oo), 
the limit uq ^ can then be taken in (10, 11). 

In order to get a continuum description of the master equation we use the expansion, 
valid when k is large, Pk+i/Pk = 1 — P/k + const. /fc^ H . We also set 

=1 + Ae + ---, (12) 



Eoo £ 
m=0 /' 



since this quantity goes to 1 for long times. The amplitude A will be determined 
below. 

For small e (i.e. large t) and fixed u (hence large fc) one gets 

4(2^ + ug') = g" + ^g' - Ag' - ^ + 0(e), (13) 
s u U 

which is in separable form. The solution for e is const. (i — fo)~^^^ or, after a change 
in the origin of time, £ = 1/ Vi. In this regime one therefore obtains 

/+(i„-^ + 2)»'+(l-|)» = «, (14) 

This equation is singular at u = and m = oo. At 7i = one expects a power 
singularity for g of the form . The Probenius series g(u) = UnU^'^'^ carried into 
(14) gives the recursion relation 

s + 71 + 2 

o„+2(s + n + l)(s + n + /3+2)-A(s + n+l)a„+i + o„ = 0, (n > -2) (15) 

with a_i = a_2 = 0. The radius of convergence of the series thus obtained is infinite. 
For n = — 2 one has (s — l)(s + /?) = 0. Hence, when u ^ 0, a basis of solutions reads 

gf\u)^u, gf\u)^u-0. (16) 

Only the first solution may be retained because of the normalisation of g (cf eqs. (10- 
11)). When u oo, an asymptotic study of the dominating terms in eq. (14) leads 
to another basis of solutions (reminiscent of the particular case A = /3 = 0, for which 
solutions are ue~" 1^ and ue~" J u~^e" /^du), namely 

gf\u) ~ ui-^e-"'/*+^«, gf\u) - (17) 

Again the second behaviour should be excluded, because of the normalisation condi- 
tions (10-11). As the two solutions must be connected, 

gf(u) = C^{A,!i)g't°\u)+G-,{A,!i)gt\u\ (18) 

and one must impose C2(^, /3) = 0. This condition determines A as a fimction of /9. 
In practice, this may be done numerically either by reconstructing g from its series or 
directly from a numerical integration of the differential equation (14). The solution 
of eq. (14) may be seen as a continuous deformation of we"" as A and /3 increase 
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Figure 2. Solution of the differential equation (14) with A = 1.9, and Monte-Carlo 
simulation (t = 1000, N = 1000, M = 500). {p = 2, (} = 4.) The numerical solution of 
(5) for large times is hardly distinguishable from the solution of (14). 

from 0. Note that the solution g{u) — we^" of the equation g" + ug' /2 + g = 0, 
corresponding to A = /3 = 0, appears also in the scaling regime of model B at T = 0, 
where fk = t~^g{kt~^^^) [3]. Figure 2 displays on the same plot the solution of eq. 
(14) with A determined following this technique, the asymptotic form of (5) for large 
times and the Monte-Carlo simulation. The agreement between all these curves is 
excellent. 

For theoretical purposes, it may be convenient to cast eq. (14) into its Schwarzian 
form 

1 / 1 . f3 (3 A^\ Ap f3il3 + 2)\ 

obtained by setting g = vw and choosing v such that no first derivative of w appears 
in the equation. This leads to 

giu)^u-^/^cxp(^-^ + ^^w{u). (20) 

Eq. (19) may be recognized as a Schrodinger equation with zero energy. Since w is 
positive, it corresponds to the ground state solution of the equation. Imposing that 
the energy of the ground state be zero determines A as a function of /3. 

Finally one may find the explicit form of the fk for small fc, i.e. in the fluid part of 
the distribution. Substituting the form (9) for fk into the master equation (5) leads 
to 

wi = "^n + ^ 

(21) 

Pk+iVk+i + PkVk-i - {Pk+i + Pk)vk = A{pk+i - Pk), {k > 0) 

the solution of which is Vk = vq + kA. The determination of vq is made possible by 
eq. (10) which gives 

poo 

vo + Apc+ / g{u)du = 0, (22) 
Jo 
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where A is already known from above, showing that vq is negative. 

Two-time correlations of the energy. At T = 0, the energy correlation function of 
model B exhibits aging [5]. It is therefore natural to expect the same property for 
model B' in the whole low temperature phase, p being fixed. Following the notations 
of [7], the correlation function of the energy of a generic box, say box number 1, at 
two times s and t {s < t) reads c{t,s) = /o(s)(.go(^i s) — /o(0)i where gk{t,s) is the 
probability that box number 1 contains k particles at time t, knowing that it was empty 
at s. The evolution in time of this quantity is given by (5), with initial conditions 
(7fc(s,s) = 5k fl. A numerical integration of these equations shows that the normalised 
correlation function c{t, s)/c, s, s) has the same asymptotic scaling form \fajt as model 
B [5], for /? > /3c, i.e. it exhibits aging in the low temperature phase. We will come 
back to the theoretical analysis of this result in a forthcoming publication. 

Discussion. We wish to emphasize the remarkable result obtained in this work, 
namely that the condensate acquires a universal scaling form directly related to the 
assumption of a regular power law behaviour of the at infinity. We checked that 
the master equation with Metropolis rule lead to the same scaling form for g{u). As 
a consequence of this scaling form, the condensation time behaves as the squared size 
of the system. 



We wish to thank R. Balian, R. Conte, J.-M. Luck and G. Mahoux for interesting 
conversations. 
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